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Following the 12 GeV upgrade, a dedicated experiment is planned with the Hall B CLAS12 detector at Jefferson 
Lab, with the aim to study electroproduction of nucleon resonances at high photon virtualities up to = 
12 GeV^. In this work we present a QCD-based approach to the theoretical interpretation of these upcoming 
results in the framework of light-cone sum rules that combine perturbative calculations with dispersion relations 
and duality. The form factors are thus expressed in terms of A* (1535) light-front wave functions at small 
transverse separations, called distribution amplitudes. The distribution amplitudes can therefore be determined 
from the comparison with the experimental data on form factors and compared to the results of lattice QCD 
simulations. The results of the corresponding next-to-leading order calculation are presented and compared 
with the existing data. We find that the form factors are dominated by the twist-four distribution amplitudes that 
are related to the P-wave three-quark wave functions of the A* (1535), i.e. to contributions of orbital angular 
momentum. 
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I. INTRODUCTION 

It is generally accepted that studies of baryon form factors 
at large momentum transfer give access to the light-front 
wave functions at small transverse separations between the 
constituents, called hadron distribution armlitudes (DAs), al¬ 
though perturbative QCD factorization (ll-l^ does not seem to 
be applicable for realistic accessible in current or planned 
experiments. The problem is that the leading contribution in¬ 
volves two hard gluon exchanges and is suppressed by the 
small factor (cts/Tr)^ ~ 0.01 compared to the “soft” (end¬ 
point) contributions which are subleading in the power count¬ 
ing in 1/Q^ but do not involve small coefficients. Hence 
the collinear factorization regime is approached very slowly. 
Model calculations suggest that “soft” contributions play the 
dominant role at present energies. Taking into account soft 
contributions is challenging because they involve a nontrivial 
overlap of nonperturbative wave functions of the initial and 
the final state hadrons, and are not factorizable, i.e. cannot be 
simplified further in terms of simpler quantities. 

In this situation the question what exactly do we learn from 
the studies of form factors is far from trivial. One exist¬ 
ing description is to introduce more complicated, transverse- 
momentum dependent (TMD) quark distributions, taking ad¬ 
vantage of Sudakov suppression of large transverse separa¬ 
tions, following the technique suggested initially by Li and 
Sterman |3l for the pion form factor. Another approach that 
we advocate in this work, is to calculate the soft contribu¬ 
tions to the form factors as an expansion in terms of nucleon 
DAs of increasing twist using dispersion relations and dualiw. 
This method is known as light-cone sum rules (LCSRs) Jj] 
and provides one with the most direct relation of the hadron 
form factors and DAs that is available at present, with no other 
nonperturbative parameters. 

The LCSR approach has been used successfully for the 
calculations of pion electromagnetic and also weak B-decay 


form factors, see Refs. |@-0] for several recent state-of-the-art 
calculations. The LCSRs for baryon form factors are more 
complicated and recent. The first application was for the nu¬ 
cleon electromagnetic form factors in Refs. Several 

further studies aimed at finding an optimal nucleon interpo¬ 
lation current CMl and extending this technique to other 
elastic or transition form factors of interest. LCSRs for the 
axial nucleon form factor were presented in lllOlHa[f^, for 
the scalar form factor in Cl and tensor form factor in ll5[] . A 
gen eralization to the full baryon octet was considered e.g. in 
d. Application of the same technique to A 7 A transitions 
was suggested in dd and to pion production at threshold 
in d- LCRSs for weak baryon decays At etc. 

were studied in iflld . etc. In the early work only the lead¬ 
ing order (LO) contributions to the coefficient functions in the 
LCSRs have been taken into account. The first complete next- 
to-leading order (NLO) analysis was done for the electromag¬ 
netic nucleon form factors in Ref. ll2l and the results appear 
to be consistent with the constraints on nucleon DAs from lat¬ 
tice calculations ll2l . The picture emerging from these stud¬ 
ies sugggests that the momentum fraction distribution of the 
valence quarks in the proton is rather broad, with ~ 40% of 
the momentum carried by the it-quark that carries proton he- 
licity, and approximately symmetric to the interchange of the 
remaining quarks. 

Our study is motivated by the dedicated experiment planned 
with the Hall B CLAS12 detector at Jefferson Lab following 
the 12 GeV upgrade, with the aim to study electroproduc¬ 
tion of nucleon resonances at high photon virtualities up to 
= 12 GeV^ jm. The corresponding form factors can 
be calculated using the LCSR machinery in terms of the DAs 
of nucleon resonances. Turning this relation around, infor¬ 
mation on the DAs of resonances can be extracted from the 
comparison of the LCSR calculations with the experimental 
data on form factors and compared to the constraints that can 
come, eventually, from lattice QCD simulations. This pro- 
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gram was suggested in Ref. ll^ and an exploratory study was 
made there for the particular case of electroproduction of the 
lowest negative parity iV*(1535) resonance. In our paper we 
elaborate on this proposal. Learning about quark distributions 
in nucleon resonances is an exciting possibility, since existing 
QCD calculations of resonance properties, e.g. on the lattice, 
rarely go beyond the mass spectrum. 

The case of TV* (1535) is special because the classifica¬ 
tion and the structure of the light-front wave functions for the 
states with opposite parity is almost identical. Hence the LC- 
SRs for the corresponding electroproduction form factors are 
very similar to the LCSRs for electromagnetic nucleon form 
factors. In particular the NLO expressions derived in Ref. ll^ 
can be overtaken with relatively minor modifications. A de¬ 
tailed analysis of these NLO LCSRs is the main goal of this 
work. 

The presentation is organized as follows. The electropro¬ 
duction form factors are introduced and the structure of the 
corresponding LCSRs is explained in Sect. 2. Sect. 3 contains 
a detailed numerical analysis and comparison with the exist¬ 
ing experimental data, and Sect. 4 is reserved for a summary 
and outlook. A large Appendix A contains a short review of 
the three-quark light-front wave functions of = ( 5 ) nu¬ 
cleon resonances, their relation to DAs, and also explains our 
conventions. Appendix B contains a simple parametrization 
of the dependence of the coefficient functions in the LC¬ 
SRs. 


II. ELECTROPRODUCTION EORM EACTORS AND 
LIGHT CONE SUM RULES 


The matrix element of the electromagnetic current j®™ be¬ 
tween spin- 1/2 states of opposite parity can be parametrized in 
terms of two independent form factors, which can be chosen 
as 


{N*{P')\3l^\N{P)) = um-{P')i^T,um{P). 


r. = 
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where q = P' — P is the momentum transfer. In what follows 
we use the standard notation = —q^. The helicity am¬ 
plitudes Ai/ 2 {Q‘^) and 5 'i/ 2 (Q^) for the electroproduction of 
Af* (1535) can be expressed in terms of the form factors ll^ : 


Q^Gi{Q‘^) + iriNimN* - rnN)G2{Q^) 


Ai /2 — eB 

pBC r 

S'1/2 = — {mj\! — mi\f)Gi{Q'^) + TnMG2{Q^) 




( 2 ) 


Here e = \/Aira is the elementary charge and B, C are kine¬ 
matic factors defined as 


The basic object of the LCSR approach to baryon form fac¬ 
tors US is the correlation function 


Tp{P,q) = i Jdxe-^‘^^N*{P')\T{jp{x)mm\0) (4) 


in which j represents the electromagnetic (or weak) probe and 
Pat is a suitable local operator with nucleon quantum num¬ 
bers. The N* resonance is explicitly represented by its state 
vector {N*{P')\, see a schematic representation in Fig.[T] The 



FIG. 1: Schematic structure of the light-cone sum rule for electro¬ 
production form factors. 

LCSR is obtained by comparing (matching) two different rep¬ 
resentations for the correlation function. On the one hand, 
when both the momentum transfer q^ = —Q^ and the mo¬ 
mentum P^ = {P' — q)^ flowing in the pjv vertex are large 
and negative, the main contribution to the integral comes from 
the light-cone region 0 and can be studied using the op¬ 
erator product expansion (OPE) of the time-ordered product 
T{j (x)77Ar(0)}. The singularity at —?► 0 of a particular con¬ 
tribution is governed by the twist of the relevant composite op¬ 
erator whose matrix element {N*\ ... |0) is related to the N* 
DA. On the other hand, one can represent the answer in form 
of the dispersion integral in P^ and define the nucleon contri¬ 
bution by the cutoff in the quark-antiquark invariant mass, the 
so-called interval of duality sq (or continuum threshold). The 
main role of the interval of duality is that it does not allow 
large momenta |fc^| > sq to flow through the Ty^r-vertex; to 
the lowest order 0 ( 0 ;°) one obtains a purely soft contribution 
to the form factor as a sum of terms ordered by twist of the 
relevant operators and hence including both the leading- and 
the higher-twist nucleon DAs. Note that the contribution of 
higher-twist DAs is suppressed by powers of the continuum 
threshold (or by powers of the Borel parameter after applying 
the usual QCD sum rule machinery), but not by powers of Q^, 
the reason being that soft contributions are not constrained to 
small transverse separations. 

The “plus” spinor projection (IA.8I) of the correlation func¬ 
tion (01 involving the “plus” component of the electromag¬ 
netic current can be parametrized in terms of two invariant 
functions 


B 


I -F (ttiat* -F toat)^ 


2m%{ 


On* 





(Q 2 


— m 


N* 


m 




4:Q‘^m 


2 

N* 


A+T+ = p+ {ruNAiQ^ P'^) -F P'^)] N+{P ), 

(5) 

where = —q^ and P'^ = {P — q^. The correlation func- 
tions A{Q‘^, P'^) and B{Q^, P'^) can be calculated in QCD 
in terms of N* DAs for sufficiently large Euclidean momenta 


C = 
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> 1 GeV^ using OPE. Schematically, 

A{Q\P'^) = fc]afe(g^P '2 Xi, 4 )Ffe(x„ 4 ), 

k 

B{Q\P'^) = Y, j[dx\bk{Q\P'^x,,^il)Fk{x,,^il)i6) 

k 


where sq — (1-5 GeV)^ is the interval of duality (also called 
continuum threshold). Matching the two above representa¬ 
tions and making the Borel transformation that eliminates sub¬ 
tractions constants 


1 


s- P'2 




(9) 


where the sum goes over all existing DAs, Fk £ 
{Vk,Ak,Tk,Sk,Pk) defined in Eq. ( |A.211 i, the integra¬ 
tion goes over quark momentum fractions and stands 
for the factorization scale. The coefficient functions 
ttkiQ^, P'^,Xi,fj.p) and bkiQ^, P '^, ) are known to the 

NLO accuracy for twist-three and twist-four DAs ll^ . and 
to leading order (LO) for twist-five and twist-six. In princi¬ 
ple this expansion also contains contributions of four-particle 
DAs with an additional gluon, five-particle with two gluons 
or a quark-antiquark pair, etc. Such contributions start at 
twist-four and they are not included in the present calcula¬ 
tion because the corresponding DAs are very poorly known 
(see, however. Ref. j^ ). It turns out that the coefficient func¬ 
tions are the same for the states with negative and positive 
parity, iV*(1535) and the nucleon, if the definitions are cho¬ 
sen as explained in Appendix A. Thus we are able to use the 
NLO expressions for the electromagnetic nucleon form fac¬ 
tors obtained in with trivial modifications, e.g. replacing 
nucleon mass mpi by m^v*. One difference is that, because 
of the larger mass, corrections of the type become 

much larger and numerically significant. Eor this reason in 
this work we use complete expressions for the LO coefficient 
functions from Ref. lUffl rather than the corresponding ex¬ 
pressions from Ref. where the expansion in powers of 
was truncated to match the accuracy of the calcu¬ 
lated NLO corrections. 

The results of the QCD calculation in Euclidean region can 
be presented in the form of a dispersion relation 

^QCD(q 2 ^ 1 

'X JQ S — P 

^QCD (q 2 p/ 2 ) ^ 1 f Im (Q^ s) + .. ( 7 ) 

Jo s — F 

where the ellipses indicate possible subtractions. The same 
correlation functions can be written in terms of physical spec¬ 
tral densities that contain a nucleon (proton) pole at P'^ —?> 
m%, nucleon resonances and the continuum. The nucleon 
contribution is, obviously, proportional to the electroproduc¬ 
tion form factors of interest, whereas for higher mass states 
one can use quark-hadron duality; 

^phys(g 2 ^p/ 2 ) 
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one obtains the sum rules 
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The dependence on the Borel parameter is unphysical and 
has to disappear in the full QCD calculation. It can be used to 
estimate theoretical uncertainties. 


III. NUMERICAL ANALYSIS 


Main nonperturbative input to the LCSRs for electropro¬ 
duction form factors is provided by the DAs of nucleon reso¬ 
nances that can be parameterized by two normalization con¬ 
stants /jv*, and a set of shape parameters (pnk, Vnk cor¬ 
responding to contributions of local operators of increasing 
dimension, see Eqs. (IA.271 i. ( IA.31b . The dependence of the 
form factors on these parameters is linear so that the results 
can conveniently be presented as 


Gi(Q') = 


Af 


N* 
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Ai 

/tv* r joo 


9T{Qn + 9i^mmo + gi\Q^)mi 


+ A {Q ) + A {Q )'-Pio + A [Q )^ii + ■ • ■ 
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and similarly 


G 2 (qA = 

, /at* 




Af* 


Af 


9T{Q^) + 9^^{Q^)mo + 9^2AQlmi 
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( 12 ) 


where the ellipses stand for the contributions of second-order 
polynomials in the leading-twist DAs (IA.27I) . terms in tf 2 o, 
^ 21 , <f 22 - The coefficient functions and 

are given by very cumbersome analytic expressions 11231] and 
depend implicitly on the masses of N* and the nucleon, the 
continuum threshold sq, Borel parameter M^, QCD cou¬ 
pling asifip) and the factorization scale fj,p. Note that, e.g., 
/™(g 2 ) includes the sum of contributions of the asymptotic 
leading-twist DA and the corresponding Wandzura-Wilczek 
terms in the higher-twist DAs, see Appendix A and Ref. 
for more details. Note also that the DA S 4 ( IA.25b correspond¬ 
ing to the iz = 2 component of the light-front three-quark 
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Method 

Af/Af 

/iv*/Ar 

‘Pio 

yjii 

¥520 

V521 

V322 

7710 

7?11 

Reference 

LCSR(l) 

0.633 

0.027 

0.36 

-0.95 

0 

0 

0 

0.00 

0.94 

this work 

LCSR (2) 

0.633 

0.027 

0.37 

-0.96 

0 

0 

0 

-0.29 

0.23 

this work 

LATTICE 

0.633(43) 

0.027(2) 

0.28(12) 

-0.86(10) 

1.7(14) 

-2.0(18) 

1.7(26) 

- 

- 

[24] 


TABLE I: Parameters of the Af*(1535) distribution amplitudes at the scale = 2 GeV^. For the lattice results (13l only statistical 
errors are shown. The set of parameters indicated as LCSR (1) corresponds to the fit to the form factors Gi(Q^) and G2(Q^) 
extracted from the measurements of helicity amplitudes in Ref. (1^ adding the errors in quadrature. The set of parameters 
indicated as LCSR (2) is obtained from the fit to helicity amplitudes including all available data at > 1.7 GeV^ (29l - l3ll . 


wave function does not contribute to the LCSRs for our choice 
of the nucleon interpolating current. 

Calculations in this work are done for the "standard" choice 
of the specific LCSR parameters; continuum threshold sq = 
(1.5 GeV)^, Borel parameter = 2 GeV^ and factorization 
(and renormalization) scale = 2 GeV^. The dependence 
on these parameters is rather mild; in particular varying the 
Borel parameter in the range 1.5 — 2 GeV^ induces an overall 
variation of form factor of the order of 10% so that, e.g., the 
ratio G 2 /Gi is largely unchanged. 

The resonance mass corrections enter the LCSRs in a com¬ 
plicated way, as terms in and The lat¬ 

ter ones do not decrease at large momentum transfers and in 
ideal case have to be resummed to all orders. The correspond¬ 
ing expression exists for the LO LCSRs ifl^ but not for the 
NLO corrections. In order to minimize this mismatch we have 
rescaled the 0{as) contributions calculated in by the ra¬ 
tio of the corresponding LO terms calculated with account for 
corrections and putting to zero. For the numeri¬ 
cally important contributions this rescaling corresponds to a 
reduction of the NLO correction by 10 — 20%. 

Existing information on the DAs of negative parity reso¬ 
nances is very scarce. The results of the recent lattice cal¬ 
culation ll24l] are presented in Table 1. The most interesting 
feature of these results is that the corrections to the asymp¬ 
totic leading twist DAs have alternating signs for the lattice 
states with increasing mass. In particular the twist-three DA 
of A^*(1535) has a very small value at the origin and is ap¬ 
proximately antisymmetric with respect to the exchange of 
the two valence quarks forming a scalar “diquark”, whereas 
the DA of A^*(1535) is symmetric and similar in shape to the 
nucleon DA, see Fig. 9 in Ref. j^ . These results are still 
exploratory and have to be taken with caution because iden¬ 
tification of lattice states with particular physical resonances 
is not obvious and requires further study. Even with this un¬ 
certainty, the lattice values are very helpful as knowing the 
order of magnitude of the parameters allows one to establish 
a hierarchy of different contributions to the LCSR. 

As an illustration, the NLO LCSR result for the form factors 
at = 2 GeV^ normalized to the dipole formula 


D{Q^) 


1 

(1 + Q7a)2 ’ 


a = 0.71 GeV^ (13) 


can be written as follows: 


GNLO(q2) 

D{Q^) 


AT* r 


Af 


0.666 — 2.187710 + 0.867711 


- 0.69/Ar. - 1.76/Ar.yjio + 1.05/Ar.yjii 
-F 1.3/Ar*(F>20 + 0.66/Ar*(F>21 - 0.06/Ar*(F>22 


GNLO(q2) 

Dm 


Af 


— 0.466 -F 1.847710 + O.O67711 


— 0.82/iv. — 1.06/Ar.</3io — l.OS/Aftpii 

+ 2.6/Ar*;F’2o + l-5/Ar*y72i + 0.39/Ar*y522 


where we use a notation /at* for the ratio of twist-three and 
twist-four couplings 


In* 



0.027(2) 


12^. 


(14) 


For comparison, the similar decomposition of the form factors 
for the LO LCSRs ll^ for the same value = 2 GeV^ reads 


Dm 


Dm 


_xr 

Af 


0.816 - 2.027710 + 0.887711 


- 0.59/iv - L60 /iv.<f>io + L19/iv.<F>ii 
+ l-26/Ar*<F’20 + 0.70/Ar*<F’21 + 0.12/Ar.(^22 


_Af* 

Af 


— 0.466 -F 1.847710 -F O.O67711 


— I.IO/at* — 0.78/Ar*(pio + 3.82/iv.(pii 
+ 2.9/Ar*y72o + L6/Ar*y72i + 0.28/Ar*y722 


so that the NLO corrections are significant. 

For convenience we provide a simple parametrization for 
the coefficient functions p"! (fTTI) . (fT2l i as functions of 
in Appendix B. This parametrization was obtained for the 
region of momentum transfers 2 GeV^ < < 12 GeV^ and 

should not be used outside this interval. In particular we found 
that the mass corrections ^ become very large for 

< 2 GeV^ so that the LCSRs become unstable (and not 
reliable). In general, different contributions to the LCSRs are 
distinguished by their Q^-dependence so that one needs a suf¬ 
ficient lever arm in to determine several of them simulta¬ 
neously. 
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FIG. 2: Helicity amplitudes 2 I 12 and S 12 for electroproduction of iV*(1535) (left panel) and the form factors 
normalized to the dipole formula (right panel). Experimental data on the left panel are taken from (s^l (empty squares) Ull (filled 
squares) (filled circles) and (triangles). The form factors on the right panel are calculated from the data ( 2 ^ on helicity 
amplitudes adding the errors in quadrature. The curves show the results of the NLO LCSR fit to the form factors Gi{Q'^) and 
G 2 (Q^) for > 1.7 GeV^ with parameters of the Af*(1535) DAs specified in the first line in Table 1 . 




0 2 4 6 8 10 12 0 2 4 6 8 10 12 

Q2 


FIG. 3 : The same as in Fig.[ 2 ]but for the fit to helicity amplitudes A12, S12 including all available data at > 1.7 GeV^. The fitted 
parameters of the A^*(1535) DAs are specified in the second line in Table 1 . 


Since the existing data for > 1.5 — 2 GeV^ are very lim¬ 
ited, we put in this work all second-order coefficients in the 
leading-twist DAs to zero, (p 20 = ¥^21 = < 7^22 = 0, used cen¬ 
tral lattice values for and , and constrained (pio, T’li 
to the lattice values within the given error bars. In this way we 
are left, essentially, with two free parameters — pio and pn. 
We expect that much more data will become available after 
the 12 GeV upgrade at Jefferson Lab where a dedicated exper¬ 
iment is planned to study electroproduction of nucleon reso¬ 
nances at high photon virtualities up to = 12 GeV^ ll^ . 

Information on the electrocouplings of nucleon resonances 
at large momentum transfers is obtained by studying electro¬ 
production of TT and p mesons in the respective resonance re¬ 
gion ll29l - [3^ . The results are usually presented for the he¬ 
licity amplitudes, and in earlier work only the larger one, 
Ai 2 (Q^), was studied for large momentum transfers. The 
latest study also includes the results on Si 2 (Q^) up to 
= 4.16 GeV^ allowing us to extract from these data the 
Dirac-like and Pauli-like transition form factors Gi(Q^) and 
G 2 {Q‘^) O that are more relevant for QCD studies. In this 


extraction we assumed that the errors for helicity amplitudes 
given in Ref. ll^ are uncorrelated and added them in quadra¬ 
ture. The results are shown in Fig. |2] and Fig. |3] on the right 
panels; it is seen that the Pauli-like form factor changes sign 
and becomes negative at large Q^, although the errors are 
quite large. 

Two different LCSR fits of the experimental data are shown 
in Figs. 1313 The difference is that in Fig. |3the fit is done 
to the form factors extracted from the data on helicity ampli¬ 
tudes reported in Ref. ll^ . and in Fig. [3 we make a fit to the 
data on helicity amplitudes Ai 2 {Q^) and Si 2 {Q^) themselves 
including all existing data for > 1.7 GeV^. In the sec¬ 
ond case the fit is driven by the data ll30l - [^ on Ati2(Q^) that 
have smaller errors and not entirely consistent with so 
that a worse description of the form factors in this fit is not a 
surprise. The corresponding parameters are listed in Table 1. 

Because of the small value of the leading twist normaliza¬ 
tion constant suggested by lattice calculations (fl^ . the results 
for Ai 2 {Q'^) and Gi{Q‘^) prove to be almost insensitive to the 
leading twist DA of the W*(1535) resonance and are domi- 
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nated by the twist-four contributions corresponding to the P- 
wave parts of the three-quark light-front wave functions (see 
Appendix A). Moreover, sensitivity of the results to the shape 
parameters of the twist-four DAs, 7710 and 7711 , is rather mild, 
cf. two last columns in the first and the second line in Table 1. 
Thus Ai 2 {Q^) and Gi(Q^) are both sensitive mostly to the 
ratio of the normalization constants which we fix to 

the lattice value 0.633 ll24l] . The LCSR predictions for these 
observables are very stable, and the agreement of the existing 
data with the normalization suggested by lattice calculations 
is encouraging. 

The Si 2 {Q^) amplitude and especially the Pauli-like form 
factor G 2 {Q^) are much more sensitive to the nonperturbative 
input and in particular to the shape parameters of the twist- 
four DAs, compare Fig. |2] and Fig. |3] Also the leading-twist 
contributions play some role in this case because of strong 
cancellations. More precise data and a larger interval in 
are needed to make this comparison quantitative. 

IV. CONCLUSIONS AND OUTLOOK 

In this work we argue that the LCSR approach can provide 
one with quantitative information on the wave functions of nu¬ 
cleon resonances at short distances. The basic idea behind this 
technique is that soft Feynman contributions to the form fac¬ 
tors are calculated in terms of small transverse distance quan¬ 
tities using dispersion relations and duality. The form factors 
are thus expressed in terms of light-front wave functions at 
small transverse separations, called DAs, without additional 
parameters. Alternatively, the distribution amplitudes can be 
extracted from the comparison with the experimental data on 
form factors and compared to the results of lattice QCD sim¬ 
ulations or other nonperturbative approaches based on, e.g., 
QCD sum rules or Dyson-Schwinger equations. The results 
of the corresponding NLO calculation for the particular case 
of the N*{1535) resonance are presented and compared with 
the existing data. We find that the form factors are dominated 
by twist-four DAs that are related to the P-wave three-quark 
wave functions, i.e., to the distribution of orbital angular mo¬ 
mentum. 

Interestingly enough the LCSRs have the same form for 


spin- 1/2 resonances of both parities so that apart from the (cal¬ 
culable) effects of resonance mass corrections the difference 
in observed form factors of, say, W*(1535) and W*(1650) can 
be attributed to the difference in the wave functions, which 
is of major interest. The differences between nucleon elastic 
form factors and electroexcitation of the Roper resonance can 
be studied in a similar manner; however, it is likely that in the 
latter case interpretation of the results may require a better un¬ 
derstanding and more sophisticated models of twist-five DAs 
than are available at present. 
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Appendices 

Appendix A: Light-front wave functions and DAs of = (A) 
nucleon resonances 

In the light-front description 131 a hadron is represented by 
the superposition of Fock states with different number of par- 
tons. Restricting ourselves to the three-quark (valence) com¬ 
ponents we view, e.g., the proton with positive helicity as a 
superposition of states with different values of the quark or¬ 
bital angular momentum projection on the direction of mo¬ 
tion, iz = — 1 , 0 , 1 , 2 , 

= (A.i) 

A nonzero value of accounts for the mismatch between the 
proton helicity and the sum of helicities of the valence quarks 
\i so that 1/2 = Ai-|-A 2 + A 3 -|-f^. The four different 
contributions can be written in terms of six independent scalar 
light-front wave functions as 




gijfc 

76 

gijfc 

7 f 

gijfe 

76 

7f 


{dx][dk±] 2 , 3) -f 2, 3) |0), 

(1,2,3) + (1,2,3)] (&I;(1)772)^^!(3) - !«) > 


y/XiX2X3 . 

[dx] [dk± 

y/XiX2X3 I 

[dx] [dfc_L 

y/XiX2X3 L' 

[dx] [dfc_L 
.^XiX2X3 L' 


fci^7/7')(l,2,3)]&l;(l)(6l;(2)&t^^(3) - b%{2p^^{3))\0) , 
'k^ki^<-^\l,2,3)] {b^:^{l)bi\{2)bll{3) - bl:^{l)bl{2)bll{3))\0). 


(A.2) 


Here ^^(l) etc. are creation operators for the quarks of spe- ment (1) stands for the dependence on longitudinal momen- 
cific flavor with positive f or negative j, helicity; the argu- 
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turn fractions and transverse momenta of the given quark, i.e. 

and so on. We use the notation for 

transverse momenta 

fc"*" = + iky , k-^ = k^ — iky . (A.3) 

The light-front wave functions (1,2,3) depend on 
momentum fractions xi and transverse momenta squared 
\k±^i\’^ = kiki of all partons. The integration measure is cho¬ 
sen as 

3 

[dx] = dxk5{l - , (A.4) 


and 




The proton light-cone DAs, in turn, are defined as matrix el¬ 
ements of gauge-invariant nonlocal operators with the three 
quark fields separated by a light-like distance. Standard de¬ 
composition 11^ involves 24 invariant functions; 


4(0|e*^''u^(ain)u;^(a2n)d!;(a3n)|W(P, A)) = 

= S'f mNCafi + ^2 ruNCafi mN {l5C)^p -f P2 ruN (u^)- 

+Vi^ {j^C)^p (75U^)^ + \vi^mN (.ixC)^^ 

+ \vrmN ( 7 -LC')a /3 ( 7 -^ 75 ^^).^ + ^ 5 ^^ (f^C')a /3 ( 75 ^^).^ + ^^^ 6 ^ (l^C')a /3 ( 75 ^^)^ 


2 pn 


+A^ {u+)^ + (jll5C)^p (mjv)7 + 2 ^^(7-l75C')«;3 (7^^^).^ 

+ ^A^ruN (7±75C')„/3 + (7^75C')„^ iu+)^ + 

+T^ {ia±pC)^^ ( 7 “^ 75 M^)., + (icr^pC')^^ ( 7 -^ 75 ^^)^ + {i<^pnC)^p ( 75 ^^), 


pn 


irrii, , I ... ml 


(*<^npC')a/3 (75WAr)^ + ^5^^ {i^XnC)^p ^ (7-^75WAr) 

+ ]-mNT^ {<T±±'C)^p + ^mArTg^ (cr±±'C')„^ (ct"^, 


7 


(A.6) 


r 


In this expression a, /3 ,7 are spinor indices, is an auxiliary 
light-like vector, = 0 , 

1 ml, r. 

= (A.7) 

where P^ is the proton momentum, P^ = m^. Further, 
= A±MAr(P, A) where un{P, A) is the usual Dirac spinor 
in relativistic normalization, the projectors are defined as 


A+ = 




9 


/4I^ 


= 9pu - 


Pp.ny + p,jn^ 

pn 


(A.8) 


and C is the charge-conjugation matrix. We use a shorthand 
notation 07 -*- = cr^i,n'^p^“ 07 Q,, etc. The invariant func¬ 
tions F = Vi,Ai,Ti correspond to contributions of a given 
collinear twist and can be written as Fourier integrals 


F{aj,pn)= y’[da;]e-*(P”)^‘^^“^P(a:,) (A. 


9) 


where F{xi) depend on the three valence quark momentum 
fractions Xi. 

Using various symmetries these functions can be combined 
in 8 independent light-cone DAs ll^ . There exists a single 
DA for the leading twist-three dH] 

(0|e*^'' (^uj{ain)Cfiuj{a2n)^ ^dJ^(a 3 n)|iV(P, A)) 

= -^fN (jm)i w]y (P) /[da;] e-^'P'A E ^ ^ 

(A. 10) 


such that m 

Ui(l,2,3) = ^fN[pN{l,2,3) + pNi2,l,3)], 
Ai(l,2,3) = i/^[(^Ar(2,l,3)-(^Ar(l,2,3)] , 
ri(l,2,3) = i/^[(^Ar(l,3,2) + (^w(2,3,l)](A.ll) 
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and for twist-four there are three independent DAs ll^ 


(^ul{ain)C'^uj{a2n)^^dl{a3n)\N{P, \)) 


X 




N,WW 

4 


(Xi) + Af (Xi) 


(A.12) 


where {xi) and are the so-called 

Wandzura-Wilczek contributions that can be expressed in 
terms of the leading-twist DA ipN{xi) The constants 

In, Af^ and A^ are defined in such a way that the integrals of 
the DAs (fiN, ^ 4 , T' 4 , S 4 are normalized to unity: 


{0\e"^^ (^ul{ain)C'^j±fuj{a2n)^ j-^'^dlia3n)\N{P, X)) 
= - i {pn)i^mNul^{P)j[dx] e-*(P") 


X 




N,WW 

4 


(xi) - Af T'f (xi) 


(A. 13) 


[dx] F{xi) = 1 . 


F € {v3Ar,$4,^4,S4}. (A.15) 


(^\[ain)Cf7/iu^^{a2n)^ 'i/id\{a3n)\N{P, A)) 

\N r 

=^{pn)'^mNuli{P)J[dx]e~'‘‘^’^^'>^^<^°-' Sf (xi), 

(A. 14) 


Using the canonical expansion of the quark fields in 
(IA.10b . (IA.12b in terms of creation and annihilation operators 
and Dirac equation to eliminate “bad” quark field components 
it is easy to calculate the required matrix elements from the 
set of light-front wave functions in E qs. JA.2b . In this way 
one obtains for our normalization (cf. 11 ^ 14^ 1 


fN<PN{xi,X2,X3) 

[Af$4^ + /^$f’'^'^](x2,Xi,X3) 
[Af4/f-/wVpf^^](Xi,X2,X3) 

A^Sf (xi,X 2 ,X 3 ) 


-4\/6 J[dk±] 2, 3), 

-sVe 


[dk±] j_ 


XsniN 


.gVe [ ^^k^ . 
J X2mN 

_ 24 V 6 f 

J XiruN 


^i^'^TV-l + (1.2,3) 

fcl^^^;l+fc2"-V^^';2l(1.2,3) 


fc^(V'^-')(l,3,2)-V-^-')(l,2,3)) 


fc2"-(V'^-'H2,3,l)-V'^-^^(2,l,3 


(A. 16) 


so that DAs correspond to integrals over the light-front wave 
functions over transverse momenta, with some prefactors. 
One has to have in mind that these relations are somewhat 
schematic since transverse momentum integrals on the right- 
hand side (r.h.s.) are divergent and have to be regulated e.g. 
introducing a cutoff. In turn, the DAs are usually defined us¬ 
ing dimensional regularization and the MS-subtraction so that 
a matching coefficient can be necessary. Also the wave func¬ 
tion renormalization factors have to be added for the quark 
fields. The twist-four DAs include additional contributions 
from the four-particle Fock states with an extra gluon 
If these contributions are taken into account, the four-particle 
quark-gluon nucleon DAs have to be added as well 1^ . 

The complete set of nucleon DAs carries the full informa¬ 
tion on the nucleon structure, in the same manner as the com¬ 
plete basis of light-front wave functions. In practice, however, 
both expansions have to be truncated and the usefulness of a 
truncated version, taking into account either the first few Fock 


states or a few lowest twist contributions, may depend on the 
concrete physics application. 

The classification of the three-quark nucleon light-front 
wave functions in Fq. (I A. 2b can be overtaken for the negative 
parity isospin-1/2 resonances, e.g. Af*(1535), without modi¬ 
fication. The symmetry under parity transformation does not 
constrain the light-front wave functions but affects the relation 
between the wave functions of the states with opposite helic- 
ity in terms of the helicity-flipped quarks. The corresponding 
expressions can be worked out using the Jacobi-Wick trans¬ 
formation M 

Y|iV, A) = pAr(-l)i/2-^|7V, -A) (A.17) 

where Y is the parity transformation followed by a 180° ro¬ 
tation along the y-axis, and is internal parity, = 1 for 

the nucleon and t/n = —1 for A^*(1535). Thus k^ k^, 

|A^, t) , whereas for the quark states q = u,d 
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^t^lO) h3> 6^^|0), but 3> - 6 ^^ 10). Applying this transformation to the both sides of Eg. (Ia!^ one obtains 


J 


WiYuVf = -m 


^ijk r 

'veJ 

\mtT^ = vN'— 
m)tf = vN 


[dx] [dkj 


^XiX2X3 I 

[dx] [dk± 


^wp(l>2,3) bl\{l)[bl^^{2)blll{3) - blj’^{2)blli{3)'^\0) 


VQ j y/XlX2X3 . 
r [dx][dkx 


Wi)uud = 


VQ j \/X\X2X3 - 
r \dx\\dk^ 


■/ 


Ve J ^JX\X2X3 . 


kti;^j^l{l,2,3) + k^i;%{l,2,3)] (bi\{l)b%{2)b],\{3) - b],\{l)b%{2)bl!^^{3))\0) 
fci^fc^V’w^(l:2,3) (bl^^{l)bl^^{2)bl\{3) - bl^^{l)by^{2)bl^^{3)^\0). 


(A.18) 


so that, e.g., for the £z = 0 states 




(0) 


Ni 


Nt 


(A. 19) 


A Lorentz-covariant definition of the DAs of negative parity 
resonances involves some freedom. It is convenient to choose 
the definition is such a way that the coefficient functions in 


the OPE of currents dUi are the same for states of both par¬ 
ities, and also the relations between different DAs imposed 
by QCD equations of motion remain the same. As noticed in 
Ref. (lH^), this can be achieved using invariant decomposi¬ 
tion of the A^* (1535) matrix element in terms of the 75 -rotated 
quark fields 


4(75)aa'(75)/3/3'(75)77'(0|e*'^'=<,(ain)M^,(a2n)d^,(a3n)|Af*(P, A)) = mN^Cap (jbUn-)^ + ■ ■ ■ (A.20) 

where the expression on the right hand side is the same as in Eq. ( IA. 6 b with obvious replacements rriM ^ mjv* etc. Projecting 
out the 75 matrices we obtain 


4(0|e*'^''uJ,(ain)-«^(a2n)(i!;(a3n)|A7(P, A)) = 

= mN-Ca0 + S 2 ruN-Cap + Pi^ rnN- {lbC)^p + P 2 wat* i'l5C)^p (75'«)^.) 

ijiC)^p (m+.)^ - 1 / 2 ^* (un.)^ + (7±C)„^ (7-^<.)7 

+ lvf'mN^ (7^Ww.) 7 - (''^*)7 - (^^)c/3 (^^•)7 


2pn 


-Af (^75 C')c/ 3 (75W)1^.)7 - Af (^75C)„^ (75UAr.)7 + 2^3 ^ 

*mA^. (7-L75C')a;3 b^l3UN*)^ - (^T'5C')a/9 (75U)^.)7 “ ¥-l3C)^p (75Mw.)7 

-Tf - Tf + ^3^*^ {^^pnC)^p {u+,)^ 


+T, 


jsf* '^N* 


pn 


(.i^npC)^p 


W* 


2 ^^ {i<^XnC)^^ il^u+,)^ - * {i<^XnC)^p 


+ -mN-Tl^ (cr_L_L-C')„^ ( 


a^^'u+ 


N-)-y+ 2^N-Ti^ {fTx±’C)^^{a 




(A.21) 


This expression replaces the decomposition (IA. 6 b for the nu¬ 
cleon. Note that there are some minus signs and in particu¬ 
lar all three leading twist DAs Vi, Ai and Ti are defined in 


our convention with a different sign as compared to the nu¬ 
cleon. As a consequence in the definition of leading-twist DA 
in terms of the chiral quark fields there is a minus sign as com- 
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pared to (lA.lOt . 

{0\e^^'^(ul{ain)Cfiu^{a2n)'j fidl{a3n)\N*{P)) 

= {pn)^ uIj. (P) f[dx] ^ (ppf, (xi), 

(A.22) 


(^uj{ain)C'^^±puj{a2n)^ ^-^'^dl{a3n)\N*{P)) 
= _ i ^pri) fi ruN-ui, (P )I[dx] ^ 


X 


/at* 'I' 


N’’,WW 

4 


(x^) - Xl'iff {x^) 


(A.24) 


where, of course, P^ = and the expressions for the 

invariant functions in terms of ipN *, are the 

same as for the nucleon, Eq. (lA.l lb . 


The twist-four DAs also acquire some signs ll^ 


(Ole*-^''^u]'(ain)C'yluj(a2n)^^d^ ( 0371)1 A^*(P)) 
= i (pn) f ul, (P )I [dx] ^ 


X 




N*,WW 

4 


(Xi)+A*$f {Xi) 


(A.23) 


( 0 |e*'^''^wJ(ain)C'^^M]'(a 2 n)^ fLdl,{a3n)\N*{p)) 

(P) j[dx] ^ Sf * (xi), 

(A.25) 


, »Af* WW / \ , T N* WW / \ • , , 

where $4 ’ (xi) and W 4 ’ (xi) are given by the same 

expressions in terms of the expansion of the leading-twist DA 

Pn* (xi) as for the nucleon. 


I 

The price to pay for universality of correlation functions for positive and negative parities is that the relations between DAs 
and light-front wave functions in this convention acquire some signs as well. 


fN*PN*{xi,X2,X3) 

[Xf‘^f+fN<^r’'^'^]iX2,X3,X3) 
[Xf'i'f -fN'Sf’'^'^]{x3.X2,X3) 
Af*5f*(xi,X2,X3) 


= +4V& y[dfcj_]'0^il.4(l,2, 3) 


= -bsVe / 

'*■^1 k* . h 


XsmN* 1 

= -sVe / 

k* . ) 

J 

X2mN* L 

= -24^6 

r ["'•■J 7 ^ 


XiUlN* 


n 


C^V’w*; 


*1 V'iv*;! + ^2 V’y;2 (1:2,3) 




fci 


N* 


(1, 3 , 2)(1,2, 3 


(- 1 ) 
N 


;'^(2,3,1)-V'y^(2,l,3 


(A.26) 


that have to be taken into account for the interpretation of the 
results. 

Parametrization of the DAs of the resonances can be over¬ 
taken from that for the nucleon. The leading-twist DA 
PN*{xi,p) can be expanded in the set of orthogonal poly¬ 
nomials Vnk[Xi) 


00 n 

PNiXi,p) = 120X1X2X3 EE Pnk{p)Pnk{Xi) , 

n=0 k=0 


j [dx] X\X 2 X 3 PYi]i{XiyPn’k' ixi] tX , (A.27) 


such that the coefficients are renormalized multiplicatively to 
one-loop accuracy. 


fN*{p) 

Pnk{p) 


fN*{Po) 

Pnk{pi)) 


( as[p) 
\<XsiPo) 
/ asjp) 

\a4po) 


2/(3/3o) 

Tnfc /00 


(A.28) 


Here /?o = 11 — §71/ is the first coefficient of the QCD , 8 - 
function and 7 jifc are the anomalous dimensions. The dou¬ 
ble sum in Eq. (IA.27b goes over a complete set of orthogonal 
polynomials Pnk{xi), fc = 0 ,..., ti, of degree n: 


Poo = 1 1 

Pio = 21(xi-X 3), Pn = 7 (xi - 2x2 + X3), 

P 20 = ^[3(xi - X 3 )^ - 3x2(xi -I-X 3 )-I- 2 x 2 ], 

P 21 = Y(a:i - 3 x 2 + a; 3 )(a;i - 0 : 3 ), 

9 

P 22 = -[x 4 -|- 9 x 2 (xi-|-X 3 )- 12 xiX 3 - 6 x 2 +X 3 ] (A.29) 
5 

etc., and the corresponding anomalous dimensions are 


Too =0, 
32 

720 =y, 


20 8 

Tio — , Til — 2 : 

40 14 

721 = -^, 722 = -^- 


(A.30) 
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The normalization condition (IA.151 I implies that <poo = 1- 
In the main text we refer to the coefficients ipnkit^o) with 
n = 1,2,..as shape parameters. The set of these coeffi¬ 
cients together with the normalization constant /Ar(p,o) at a 
reference scale po specifies the momentum fraction distribu¬ 
tion of valence quarks on the nucleon. They are related to 
matrix elements of local gauge-invariant three- quar k opera¬ 
tors and can be calculated, e.g., on the lattice 124]^. 

The twist-four DAs can be parameterized as 11^ 

= 24x1X2^1+ r|loi^J-)^^loix3,Xl,X2) 

- riii{^)TZii{x3, xi, X2)'^, 

T'f {x^,fi) = 24a;ia;3|l-I-77io(M)^io(a;2,a;3,a;i) 

-I- 77ii(p)7^ii(a;2,a;3,a;i)|, 

Sf {Xi,n) = 24a:2a:3|l-f ^^lo(/r)7^lo(a;l,X3,a;2)|, 

(A.31) 

where 


^io(a;i,a;2,X3) = 4 \^xi + X 2 --X 3 

TZii{xi,X 2 ,X 3 ) = [xi- X 2 + \x 3 ] (A.32) 


and 7710 (/i), pii(p), ^10 (fi) are the new shape parameters. The 
corresponding one-loop anomalous dimensions are S 


in) _ ^ (n) _ 4 

7l0 g ’ ill ^ ’ 


7 ^ = y ■ (A-33) 


For the twist-five DAs we take into account contributions of 
geometric twist-three and twist-four operators as explained in 
Ref. m. 

Note that the asymptotic DAs (at very large scales) for the 
nucleon and the resonances are the same: 


'4}'^{xi) = 120x1X2X3 , 4 > 4 ®(a:i) = 24 xiX 2 , 

= 24x1X2(1 + ^(1 “ 5x3)), 

= 24x1X3(1 -I- ^(1 - 5x2)), 

S 4 (xi) = 24 X 2 X 3 , 4 'f(xj) = 24 xiX 3 . (A. 34 ) 


For completeness we also give here the definitions of the 
normalization constants in terms of matrix elements of local 
three-quark operators: 


(01 {u,Ciiu,m^3idk (0) I N\P)) 

= fN»{pn)-^3iUN-'{P), 
(0|e*'^''(uy7^Uj)(0)757''4(0)|iV*(p)) 

= Af rriN.'y^UN^iP), 
(0|e*^'=(uya^,U4)(0)75a'^"4(0)|lV*(P)) 

= X 2 mN-l3UN-{P)- (A.35) 


Appendix B: Parametrization of coefficient functions 


For convenience we provide a simple parametrization for 
the coefficient functions f^\, gi \ appearing in (fTTli . (fT2li . for 
the range 2 < < 12 GeV^: 



(B.l) 


where D{Q^) is the dipole form factor ( [13] ). The coefficients 
and are collected in Table |II] 


[1] V. L. Chernyak and A. R. Zhitnitsky, JETP Lett. 25, 510 (1977) 
Sov. J. Nucl. Phys. 31, 544 (1980); V. L. Chernyak, A. R. Zhit¬ 
nitsky and V. G. Serbo, JETP Lett. 26, 594 (1977) Sov. J. Nucl. 
Phys. 31, 552 (1980). 

[2] A. V. Radyushkin, JINR report R2-10717 (1977), 

arXiv:hep-ph/0410276 (English translation); 

A. V. Efremov and A. V. Radyushkin, Theor. Math. Phys. 42, 
97 (1980) Phys. Lett. B 94, 245 (1980). 

[3] G. P. Lepage and S. J. Brodsky, Phys. Lett. B 87, 359 (1979); 


Phys. Rev. D 22, 2157 (1980). 

[4] H. -n. Li and G. E. Sterman, Nucl. Phys. B 381, 129 (1992). 

[5] 1. 1. Balitsky, V. M. Braun and A. V. Kolesnichenko, Sov. J. 
Nucl. Phys. 44, 1028 (1986); Nucl. Phys. B 312, 509 (1989); 
V. L. Chernyak and I. R. Zhitnitsky, Nucl. Phys. B 345, 137 
(1990). 

[6] V. M. Braun, A. Khodjamirian and M. Maul, Phys. Rev. D 61, 
073004 (2000); 

J. Bijnens and A. Khodjamirian, Eur. Phys. J. C 26, 67 (2002). 







12 




<1 


lOO 


^p;l 

l20 
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l22 

o 

II 
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p = i 

0.773867 

-6.0864 

0.646869 

-0.307557 

-1.94975 

-0.706607 

4.61356 

3.7677 

0.112153 

p = 2 

-0.18913 

5.62993 
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